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ABSTRACT 


Two  complementary  methodologies  are  described  to  quantify  the  effects  of  crack-tip  stress 
triaxiality  {constraint)  on  the  macroscopic  measvires  of  elastic-plastic  fracture  toughness, 
J  and  Crack-Tip  Opening  Displacement  (CTOD).  In  the  continuum  mechanics  methodolo¬ 
gy,  two  parameters,  c/and  Q,  suffice  to  characterize  the  full  range  of  near-tip  environments 
at  the  onset  of  fracture.  J  sets  the  size  scale  of  the  zone  of  high  stresses  and  large  deforma¬ 
tions  while  Q  scales  the  near-tip  stress  level  relative  to  a  high  triaxiality  reference  stress 
state.  Fiill-field  finite  element  calculations  show  that  the  J-Q  field  dominates  over  physi¬ 
cally  significant  size  scales,  i.e.,  it  describes  the  environment  in  which  brittle  and  ductile 
failure  mechanisms  are  active.  The  material’s  fracture  resistance  is  characterized  by  a 
toughness  locus,  t7c(Q).  which  defines  the  sequence  of  J-Q  values  at  fracture  determined 
by  experiment  firom  high  constraint  conditions  (Q  =®0)  to  low  constraint  conditions  (Q  <  0). 

To  reduce  experimental  effort  needed  to  construct  a  J-Q  toughness  locus,  a  micromechanics 
methodology  is  described  which  predicts  the  toughness  locus  using  crack-tip  stress  fields 
and  critical  J-values  from  a  few  fracture  toughness  tests.  A  robust  micromechanics  model 
for  cleavage  fi:^cture  has  evolved  from  the  observations  of  a  strong,  spatial  self-similarity 
of  crack-tip  principal  stresses  under  increased  loading  and  across  different  fracture  speci¬ 
mens.  While  the  spatial  variation  remains  self-similar,  the  magnitudes  of  principal 
stresses  vary  dramatically  as  crack-tip  constraint  evolves  under  loading.  The  microme¬ 
chanics  model  employs  the  volume  of  material  bounded  within  principal  stress  contours  at 
fracture  to  correlate  Jc  values  for  different  specimens  and  loading  modes.  The  J-Q  descrip¬ 
tion  of  the  crack-tip  stress  fields  predicts  the  similarity  of  principal  stress  contours  as 
constraint  evolves  under  loading.  For  an  applied  J-value,  the  size,  but  not  the  shape,  of 
principal  stress  contoxirs  is  altered  by  the  near-tip,  uniform  hydrostatic  stress  states  of  ad¬ 
justable  magnitude  characterized  by  Q.  These  observations  imply  that  values  specified  for 
metallvu-gical  parameters  in  the  micromechanics  model,  such  as  the  critical  fracture  stress 
and  the  distance  to  the  critical  particle,  have  only  a  weak  influence  on  the  relative  variation 
of  fracture  toughness,  Jc,  with  constraint  for  a  given  material  and  temperature. 

This  report  explores  the  fundamental  concepts  of  the  J-Q  description  of  crack-tip  fields, 
the  fracture  toughness  locus  and  micromechanics  approaches  to  predict  the  variability  of 
macroscopic  fractxxre  toughness  with  constraint  under  elastic-plastic  conditions.  While 
these  concepts  derived  from  plane-strain  considerations,  initial  applications  in  fully  3-D 
geometries  are  very  promising.  Computational  results  are  presented  for  a  siirface  cracked 
plate  containing  a  6:1  semi-ellipticsd,  a=f/4  flaw  subjected  to  remote  uniaxial  and  biaxial 
tension.  Crack-tip  stress  fields  consistent  with  the  J-Q  theory  are  demonstrated  to  exist 
at  each  location  along  the  crack  front.  The  micromechanics  model  employs  the  J-Q  descrip¬ 
tion  of  crack-frant  stresses  to  interpret  fracture  toughness  values  measiired  on  laboratory 
specimens  for  firacture  assessment  of  the  surface  cracked  plate.  The  computational  results 
suggest  only  a  minor  effect  of  the  biaxial  loading  on  the  crack  tip  stress  fields  and,  conse¬ 
quently,  on  the  propensity  for  fracture  relative  to  the  uniaxial  loading. 
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Continuum  and  Micromechanics  Treatment  of 
Constraint  in  Fracture 


1.  INTRODUCTION 

TVo  fundamental  concepts  underlie  both  linear-elastic  fracture  mechanics  and  elastic- 
plastic  fracture  mechanics:  [1]  the  relevant  crack-tip  singularity  dominates  over  micro- 
structurally  significant  size  scales  and  [2]  the  parameter  iiT/,  or  J,  uniquely  scales  the  ampli¬ 
tude  of  the  near-tip  fields.  In  an  actual  structure,  the  crack-tip  field  must  be  perturbed  by 
the  external  boundary  and  the  loading  distribution  {from  afar),  and  by  the  zone  of  inelastic¬ 
ity  and  small-scale  heterogeneities,  e.g.,  grains,  microcracks,  crack  face  roughness  {from 
within).  However,  when  the  zone  of  inelasticity  and  small-scale  heterogeneities  remains 
small  compared  to  the  external  geometry,  the  asymptotic  field  is  approximately  unper¬ 
turbed  in  an  annulus  which  is  larger  than  the  zone  of  inelasticity  and  small-scale  heteroge¬ 
neities,  but  much  smaller  than  the  external  geometry.  Strain-stress  fields  in  such  an  annu¬ 
lus  are  determined  completely  by  the  singularity  solution  (see  review  article  by  Hutchinson 
[1]).  Under  this  condition,  the  effects  of  remote  loading  and  external  boundaries  are  com¬ 
municated  to  the  crack  tip  through  Kj,  or  *7,  alone.  Moreover,  when  linear  elasticity  prevails 
at  the  macro-scale,  then  boundaiy  loading  as  well  as  traction  free  boundaries  exert  their 
influence  on  the  near-tip  field  through  Ki  alone,  with  no  effect  on  the  actual  distribution. 
Similarly,  the  plastic  fields  for  well-contained  yielding  sense  external  boundaries  and  load¬ 
ing  only  through  J  [2,3,4].  However,  for  large-scale  yielding  in  finite  bodies,  ♦  n- 

ship  between  the  scaling  parameter,  J,  and  the  near-tip  fields  loses  the  one-to-oUw  corre¬ 
spondence  [5,6,7]  This  loss  of  uniqueness,  often  termed  loss  of  constraint,  produces  the  in¬ 
creases  in  fi:ncture  toughness  observed  for  tension  geometries  and  for  shallow  notch  bend 
specimens.  The  mismatch  of  constraint  conditions  at  the  crack  tip  appeu'ently  plays  a  domi¬ 
nant  role  in  the  often  disappointing  correlation  between  fracture  specimen  b^avior  [C(T), 
SE(B)]  and  the  behavior  observed  in  large-scale,  tension  loaded  tests. 

Constraint  effects  are  most  pronounced  for  low-to-medium  strength  structural  steels 
(and  their  weldments)  operating  in  the  ductile-to— brittle  transition  region  where  unstable 
firacture  occurs  by  the  micromechanism  of  transgranular  cleavage.  Many  nuclear,  civil  and 
marine  structures  operate  in  the  transition  region  over  significant  portions  of  their  life¬ 
times.  Unlike  the  more  ductile  mechanism  of  slow  stable  tearing,  cleavage  fractures  most 
often  trigger  catastrophic  failure  of  even  highly  redvmdant  structural  systems.  Extensive 
experimental  studies  (see,  for  example  [8-12])  have  readily  demonstrated  the  much  great¬ 
er  sensitivity  of  cleavage  fracture  toughness,  Jc,  to  constraint  than  is  observed  for  ductile 
initiation  toughness,  Jjc,  and  for  ductile  crack  growth  resistance  J-Aa. 

In  the  past  three  years,  new  approaches  have  appeared  to  quantify  constraint  and  to 
predict  the  effects  of  constraint  changes  on  macroscopic  (engineering)  fracture  toughness 
characterized  by  J  and  the  crack  tip  opening  displacement,  CTOD  or  6.  Very  detailed,  elas¬ 
tic-plastic  finite  element  analyses  provide  correlations  of  crack-tip  stress  fields  over  dis¬ 
tances  r  s  l-8<5  with  loading  level  (J),  loading  mode  (tension  vs.  bending),  specimen  geome¬ 
try  and  strain  hardening.  Such  computations  stimulated  development  of  the  J~Q  continu¬ 
um  mechanics  framework  [13-19]  to  describe  the  near-tip  fields  under  very  general  condi¬ 
tions  of  loading  in  finite  bodies.  ^Afithin  this  framework,  the  J-integral  sets  the  scale  of  de¬ 
formation  at  the  crack  tip  (i.e.,  the  CTOD)  while  the  hydrostatic  stress  parameter,  Q,  quan- 
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tifies  the  level  of  stress  triaxiality  over  distances  r  <  1-86  ahead  of  the  tip  in  which  the  mi¬ 
croseparation  processes  occur.  Under  increased  loading,  each  fracture  specimen  for  a  specif¬ 
ic  material/temperature  follows  a  characteristic  J-Q  driving  force  curv’e  which  defines  the 
evolution  of  crack-tip  deformation  and  constraint.  Specimens  fracture  at  critical  J-values 
which  depend  on  Q  (critical  J-values  are  determined  by  laboratory  testing,  the  correspond¬ 
ing  Q-values  at  fracture  are  determined  from  analysis).  By  testing  fracture  specimens  that 
exhibit  a  wide-range  of  constraint  conditions  (e.g.  shallow-to-deep  notch  SE(B)  speci¬ 
mens),  the  toughness  locus  for  the  material  is  constructed,  i.e.,  the  ciirve  connecting  all  crit¬ 
ical  J  vs.  Q  points. 

The  level  of  stress  triaxiality,  quantified  conveniently  by  Q,  plays  a  dominant  role  in  the 
competition  between  fracture  initiation  by  cleavage  and  ductile  tearing  for  ferritic  materi¬ 
als  in  the  ductile-to-brittle  transition  region.  Cleavage  fracture  is  controlled  by  critical  lev¬ 
els  of  the  hoop  stress  (opening)  acting  over  microstructuraUy  significant  distances  ahead 
of  the  crack  tip  [20,21].  Ductile  tearing  is  controlled  by  several  competing  processes,  includ¬ 
ing  growth  of  a  dominant  void  and  coalesescence  with  the  crack  tip,  interaction  of  many 
growing  voids  leading  to  localization  of  plastic  flow,  and  zigzag  mode  of  ductile  tearing 
[22,23].  Each  of  these  mechanisms  is  dependent  to  some  extent  on  the  mean  stress.  Since 
Q  quantifies  both  the  hoop  stress  and  the  mean  stress  over  the  operative  length  scales  for 
both  cleavage  and  ductile  micromechanisms,  it  provides  a  common  parameter  to  interpret 
fractvu'e  initiation  in  terms  of  a  single  toughness  locus. 

Another  two-parameter  approach  also  receiving  considerable  attention  utilizes  J  and 
the  elastic  T-stress  [7,24-29].  These  studies  propose  to  correlate  crack-tip  stress  triaxial¬ 
ity  in  contained  and  fully-yielded  cracked  bodies  using  the  the  elastic  T-stress.  The  J-T 
and  J~Q  approaches  are  eqmvalent  under  well-contained  yielding  conditions  [14-16]. 
However,  under  fully-yielded  conditions  the  T-stress  becomes  undefined;  moreover,  T  is 
proportional  to  Kj  and  near  limit  load  Kj  (and  thus  T)  approaches  a  saturation  value  inde¬ 
pendent  of  additional  plastic  deformation.  In  contrast,  the  Q-parameter  continues  to 
evolve  over  the  entire  range  of  plastic  yielding.  Numerical  studies  have  shown  that  the  J-T 
approach  overestimates  the  actual  stress  triaxiality  for  some  geometries  and  underesti¬ 
mates  it  in  other  cases  so  that  there  is  not  a  consistent  trend  [16,30].  An  extensive  study 
of  the  limits  of  applicability  of  the  T-stress  as  a  correlator  of  near-tip  stress  triaxiality  can 
be  found  in  [30].  Readers  are  referred  to  the  publications  [7,24-29]  and  references  therein 
for  details  of  the  J-T  approach.  This  paper  focuses  on  continuum  and  micromechanics  ap¬ 
proaches  believed  to  have  broader  applicability. 

The  J-Q  approach  may  become  prohibitively  expensive  as  the  number  of  specimens  and 
temperatures  of  interest  increases,  lb  reduce  the  cost,  a  micromechanics  model  for  cleavage 
fracture  is  introduced  to  predict  the  toughness  locus  using  the  finite  element  stress  fields 
and  the  Jc  values  from  a  few  fracture  toughness  tests.  Recent  developments  [31,32]  in  the 
formulation  of  a  robust  micromechanics  model  focus  on  the  observation  of  a  strong,  spatial 
self-similarity  of  crack-tip  principal  stresses  under  increased  loading  and  across  different 
fractxire  specimens.  While  the  spatial  vauiation  remains  self-similar,  the  magnitudes  of 
principal  stresses  vary  dramatically  as  crack-tip  constraint  evolves  under  loading.  The  mi¬ 
cromechanics  model  employs  the  volume  of  material  bounded  within  principal  stress  con¬ 
tours  at  fracture  to  correlate  Jc  vailues  for  different  specimens  and  loading  modes.  The  simi¬ 
larity  of  principal  stress  contours  as  constraint  evolves  under  loading  is  entirely  consistent 
with  the  J-Q  description  of  the  crack-tip  stress  fields.  For  an  applied  J-value,  the  size,  but 
not  the  shape,  of  principal  stress  contours  is  altered  by  the  near-tip,  uniform  hydrostatic 
stress  states  of  adjustable  magnitude  characterized  by  Q.  These  observations  imply  that 
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values  specified  for  metallurgical  parameters  in  the  micromechanics  model,  such  as  the 
critical  fracture  stress  and  the  distance  to  the  critical  particle,  have  only  a  weak  influence 
on  the  relative  variation  of  fracture  toughness,  Jc,  with  constraint  for  a  given  material  and 
temperature. 

This  paper  explores  the  fundamental  concepts  of  the  J-Q  description  of  crack-tip  fields, 
the  fracture  toughness  locus  and  micromechanics  approaches  to  predict  the  variability  of 
macroscopic  fracture  toughness  with  constraint  under  elastic-plastic  conditions.  While 
these  concepts  derived  from  plzine-strain  considerations,  initial  applications  in  fully  3-D 
geometries  are  very  promising.  Computational  results  are  presented  for  a  surface  cracked 
plate  containing  a  6:1  semi-elliptical,  a=f/4  flaw  subjected  to  remote  uniaxial  and  biaxial 
tension.  Crack-tip  stress  fields  consistent  with  the  J-Q  theory  are  demonstrated  to  exist 
at  each  location  along  the  crack  front.  The  micromechanics  model  employs  the  J-Q  descrip¬ 
tion  of  crack-front  stresses  to  interpret  fracture  toughness  values  measured  on  laboratory 
specimens  for  fracture  assessment  of  the  surface  cracked  plate.  The  computational  results 
suggest  only  a  minor  effect  of  the  biaxial  loading  on  the  crack  tip  stress  fields  and,  conse¬ 
quently,  on  the  propensity  for  fracture  relative  to  the  uniaxial  loading. 
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2.  J-Q  THEORY 

Consider  a  cracked  bod  >  of  characteristic  dimension  L  loaded  remotely  by  a  stress  denoted 
aoo .  The  scale  of  crack-tip  deformation  is  measured  by  J/ao  where  oq  is  the  material’s  ten¬ 
sile  yield  stress  (d  «  J/ao).  At  a  sufficiently  low  load,  L>  J/ao  and  it  can  be  shown  from  di¬ 
mensional  grounds  that  all  near-tip  fields  are  members  of  a  single  family  of  crack-tip 
fields.  Each  member  field  is  characterized  by  its  level  of  deformation  as  measured  by  J /ao 
and  by  its  level  of  crack  tip  stress  triaxiedity  as  measured  by  Q,  which  also  identifies  that 
field  as  a  particular  member  of  the  family.  For  example,  the  self-similar  solution  of  Rice 
and  Johnson  [33]  and  McMeeking  [34]  (as  well  as  the  HRR  field  [3,4])  is  the  Q  s  0  member 
field.  The  Q-family  of  fields  provides  the  proper  characterizing  parameter  for  the  full  range 
of  near-tip  stress  states. 

In  the  following  discussion,  attention  is  directed  to  the  prospective  fracture  region 
ahead  of  the  crack  tip  on  the  scale  of  several  crack  opening  displacements,  6,  representing 
the  environment  in  which  the  failiire  mechanisms  are  active. 


2.1  Q-Family  of  Fields-MBL  Formulation 


The  Q-family  of  fields  is  constructed  using  a  modified  boundary  layer  (MBL)  formulation 
in  which  the  remote  tractions  are  given  by  the  first  two  terms  of  the  small-displacement- 
gradient  linear  elastic  solution  (Williams  [35]), 


(2.1) 


Here  r  and  6  are  polar  coordinates  centered  at  the  crack  tip  with  0=  0  corresponding  to  a  line 
ahead  of  the  crack  is  shown  in  the  insert  in  Fig.  2.1.  Cartesian  coordinates,  X  and  Y  with 
the  .X— axis  running  directly  ahead  of  the  crack,  are  used  when  it  is  convenient.  Within  the 
MBL  formulation, 


J  =  (2.2) 

under  plane  strain  conditions,  where  J  is  Rice’s  J-integral  [1],  E  is  Yoxmg’s  modulus  and 
V  is  Poisson’s  ratio. 


Fields  of  different  crack  tip  stress  triaxialities  can  be  induced  by  applying  different  com¬ 
binations  of  K  and  T.  From  dimensional  considerations,  these  fields  can  be  organized  into 
a  family  of  crack  tip  fields  parameterized  by  T/oq: 

=  oro4^^,e;r/aoj  .  (2.3) 

That  is,  the  load  parameter  T /cq  provides  a  convenient  means  to  investigate  and  parame¬ 
terize  specimen  geometry  effects  on  near-tip  stress  triaxiality  under  conditions  of  well- 
contained  yielding.  Such  studies  have  been  carried  out  by  Betegon  and  Hancock  [24],  Bilby 
et  al.  [36]  and  Harlin  and  Willis  [37].  Nevertheless,  the  result  in  (2.3)  cannot  have  general 
applicability  since  the  elastic  solution  (2.1),  upon  which  the  T-stress  is  defined,  is  an 
asymptotic  condition  which  is  increasingly  violated  as  plastic  flow’  progresses  beyond  well- 
contained  yielding. 

Recognizing  the  above  limitation,  O’Dowd  and  Shih  [  13,14],  referred  to  as  OS,  identified 
members  of  the  family  of  fields  by  the  parameter  Q  which  arises  naturally  in  the  plasticity 
analysis.  OS  write: 
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(2.4) 


^ij  ^0  f  ij 


^ofij 


The  additional  dependence  of  fij,  gij  and  hi  on  dimensionless  combinations  of  material  pa¬ 
rameters  is  understood.  The  form  in  (2.4)  constitutes  a  one-parameter  family  of  self-simi¬ 
lar  solutions,  or  in  short  a  Q-family  of  solutions.  The  annular  zone  over  which  (2.4)  accu¬ 
rately  quantifies  the  actual  field  is  called  the  t7-Q  annulus. 


2.2  Difference  Field  and  Near-Tip  Stress  TriaxiaUty 


Using  the  modified  boundary  layer  formulation,  and  considering  a  piecewise,  power-law 
hardening  material,  OS  generated  the  full  range  of  small  scale  yielding,  plane  strain  solu¬ 
tions,  designated  by  (uy  )ssY-  OS  considered  the  difference  field  defined  by 

AOjj  =  (<7y)ssY  “  (^y^HRR 

where  (Uy  )hrr  is  the  HRR  field.  They  systematically  investigated  the  difference  field  within 
the  forward  sector  1 0 1  <  ^/2  of  the  annulus  J  Jao<r  <  5J  Ioq  ,  since  this  zone  encompasses 
the  microstructuraily  significant  length  scales  for  both  brittle  and  ductile  fracture[20].  Re¬ 
markably,  the  difference  field  in  the  forward  sector  displayed  minimal  dependence  on  r. 
Noting  this  behavior,  OS  expressed  the  difference  field  within  the  forward  sector  in  the  form 

AOij  =  Qa^Oife)  ,  (2.6) 

where  the  angular  functions  are  normalized  by  requiring  ay(0=O)  to  equal  unity.  More¬ 
over,  the  angular  functions  within  the  forward  sector  exhibit  these  features:  arr=‘(^ee‘^ 
constant  and  \are\<\aee  \  (see  Figs.  3, 4,  and  5  in  [13]). 

Thus  the  difference  field  within  the  sector  1 0  ]  <  7il2  andJ  Ioq  <r  <5JIoo,  correspond  ef¬ 
fectively  to  a  spatially  uniform  hydrostatic  stress  state  of  adjustable  magnitude,  i.e. 
Aaij=Qaodij.  Therefore  Q  defined  by 

Q  _  ~..y  HTO  at  0  «  0,  r  =  ZJ/a^  (2.7) 

is  a  natural  measure  of  near-tip  stress  triaxiality,  or  crack  tip  constraint,  relative  to  a  high 
triaxiality  reference  stress  state.  For  definiteness  we  have  evaluated  Q  at  r=2Jlao,  however 
we  note  that  Q  is  efff actively  independent  of  distance.  In  words,  Q  is  the  difference  between 
the  actual  hoop  stress  and  the  corresponding  HRR  stress  component  at  r=2Jloo,  the  differ¬ 
ence  being  normalized  by  oq.  The  distance  chosen  for  the  definition  of  Q  lies  just  outside  the 
finite  strain  blunting  zone.  It  is  preferable  that  Q  be  defined  at  a  distance  which  is  some 
multiple  of  the  crack  tip  opening  displacement;  the  present  definition  suffices  for  our  pur¬ 
poses. 

OS  also  considered  the  difference  field  whereby  the  reference  solution  is  the  standard 
small  scale  yielding  solution,  (ui/lss^T-o  which  is  driven  by  Kj  alone,  i.e., 

Acfij  “  (o^y)sSY  “  ^‘^y)sSY;T-0  ^2.8) 

In  this  case  the  difference  field  in  the  for  ard  sector  matches  a  spatially  uniform  hydrostat¬ 
ic  stress  state  even  more  closely.  Thus  an  alternative  definition  of  Q  is 

Q  _  -  (<^w)ssY;T-o  at  0  0,  r  =  2J/aQ  .  (2.9) 

“0 

Representative  stress  distributions  of  the  Q-family  of  fields  can  be  found  in  [13,14]. 


5 


Figure  2.1  Definition  of  the  Modified  Boundary  Layer  (MBL)  problem. 


2.3  Choice  of  Reference  Field 


The  value  of  Q  is  slightly  affected  by  the  choice  of  reference  field.  Thus  a  small  increment 
(or  decrement)  must  be  applied  to  the  Q-values  if  the  reference  field  is  changed  from 
(o&fl)HRR  to  (o&e)sSY:T=Oi  or  vice  versa.  The  reference  field  distributions  according  to  the  HRR 
singularity  and  the  small  scale  yielding  solutions  for  small  strain  and  finite  strain  are  given 
in  l^ble  2.1.  The  material’s  uniaxial  stress-strain  response  is  represented  by  an  elastic 
power-law  model  having  the  form 


^0  “  O’o/^ 


(2.10) 


r  a/E  if  n  <  ffo 

^  if  (t>(Tq 

with  values  of  E/tTo=50C,  v=0.3  adopted  in  the  computations.  Figure  2.2  shows  typical  refer¬ 
ence  fields  determined  from  the  MBL  formulation  with  r=0. 


In  practice  it  really  does  not  matter  whether  we  use  (<^)hrr  (2.7),  or  (cr(9e)ssY;T=o  (2.9), 
for  the  definition  of  Q  so  long  as  it  is  applied  consistently.  In  other  words,  the  evaluation 
and  tabulation  of  Q  solutions  for  test  specimens,  the  determination  of  the  toughness  locus 
from  test  data,  and  subsequent  applications  of  such  data  to  predict  fracture  in  structural 
components  shovdd  be  based  on  the  same  reference  field.  Nevertheless,  use  of  (2.9)  can  ex¬ 
tend  the  range  of  applicability  of  the  J-Q  approach  and  is  preferable  when  it  is  desired  to 
assess  the  spatial  extent  of  the  J-Q  annulus.-  A  parameter  which  can  ascertain  the  robust¬ 
ness  of  the  J-Q  field  is  discussed  in  the  next  Section. 

A  reference  distribution  determined  from  a  small-displacement-gradient  analysis  is 
adequate  for  most  applications.  However,  accurate  descriptions  of  fields  near  the  zone  of 
finite  strains  may  be  desirable  in  some  applicatioiis,  e.g.  computational  studies  on  the  mi¬ 
cromechanisms  of  ductile  initiation.  In  such  cases  it  is  preferable  to  calculate  (uflfl)ssY;T*o 
by  a  finite  deformation  analysis  and  to  use  (2.9)  for  the  definition  of  Q.  More  importantly 
for  practical  applications,  the  evolution  of  stress  triaxiality  in  a  finite-width,  cracked  body 
can  be  evaluated  for  an  actual  stress-strain  relation,  not  just  the  power-law  relation,  if  the 
(o&fl)sSY;T*o  reference  field  for  the  MBL  model  is  determined  with  the  same  stress-strain 
relation.  This  extends  the  applicability  of  the  approach  to  a  much  broader  range  of  material 
responses.  In  contrast,  the  reference  field  (o^Ihrr  is  defined  for  an  elastic  power-law  hard¬ 
ening  material  and  the  calculations  in  the  finite  body  also  must  employ  an  elastic  power- 
law  hardening  relation. 


2.4  Variatioii  of  Q  with  Distance 


Because  Q  scales  the  difference  field  relative  to  a  reference  stress  state,  it  provides  a  sensi¬ 
tive  measure  of  the  evolution  of  near-tip  stress  triaxiality  in  finite  width  cracked  bodies. 
It  Eilso  can  be  used  to  detect  changes  in  the  stress  triaxiality  that  deviates  from  the  pattern 
that  develops  iinder  MBL  loadings.  For  this  purpose,  we  consider  Q{r)  defined  by 


Q(r) 


at  0  =  0 


(2.11) 


where  r  »  Note  that  (os^)ssY:T»0  is  chosen  as  the  reference  field. 

The  mean  gradient  of  Q  over  1  <  r  <  5, 


«•  -  VV-5)-9g-.l)  .  (2.12) 

can  be  used  to  monitor  changes  in  the  pattern  of  the  stress  triaxiality  ahead  of  the  crack 
that  do  not  conform  to  a  spatially  uniform  hydrostatic  stress  field  of  adjustable  magnitude. 
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Q'  provides  a  measure  of  the  robustness  ofthe  J-Q  fields  in  the  application  of  interest.  For 
exEunple,  Q'=0.04  means  varies  by  less  than  0.16ao  over  the  interval  1  <  r  <  5;  that 
is,zlc%fl  is  effectively  constant  over  those  distances.  On  the  other  hand,  an  IQ'  1  much  larger 
than  0.1  implies  that  the  variation  ofAagQ  over  the  interval  1  <  r  <  5  can  be  comparable 
to  CTo-  This  is  \inacceptably  large  if  the  theory  is  employed  to  predict  cleavage  fracture  w'hich 
is  very  sensitive  to  changes  in  the  hoop  stress. 


Table  2.1:  Reference  stresses,  o^/ao,  for  MBL  problem,  T/ao=0. 


n 

r/(J/{7o) 

HRR 

Small  Strain 

Finite  Strain 

1 

5.99 

5.46 

5.95 

2 

5.04 

4.53 

4.72 

3 

3 

4.55 

4.06 

4.19 

4 

4.24 

3.76 

3.85 

5 

4.01 

3.53 

3.61 

1 

4.42 

4.83 

2 

3.90 

4.06 

5 

3 

3.63 

3.73 

4 

3.79 

3.44 

3.52 

5 

3.65 

3.29 

3.36 

1 

3.93 

3??9 

2 

3.59 

3.52 

10 

3 

3.46 

3.22 

3.33 

4 

3.38 

3.12 

3.20 

5 

3.31 

3.03 

3.11 

1 

— 

2.50 

2 

• 

2.97 

00 

3 

— 

2.77 

2.91 

4 

— 

2.74 

2.86 

5 

- 

2.71 

2.82 

2.5  Simplified  Forms  for  Engineering  Applications 

Two  simplified  representations  for  the  Q-family  of  fields  within  the  forward  sector  have 
been  proposed  by  OS.  The  first  is 

<Ty  =  (CTy)HRR  +  (2.13) 

where  d/j  is  the  Kronecker  delta.  This  form  is  consistent  with  (2.7).  The  second  form  is 

^ij  ~  ^^ipSSY;T~0  Q^cAij  (2.14) 

which  is  consistent  with  (2.9).  The  physical  interpretation  of  (2.13)  and  (2.14)  is  this:  nega¬ 
tive  (positive)  Q  values  mean  that  the  hydrostatic  stress  ahead  of  the  crack  is  reduced  (in¬ 
creased)  by  Qao  from  the  J-dominant  stress  state,  or  the  standard  small  scale  yielding 
stress  state.  This  interpretation  is  precise  when  |  Q'  |  <<  1. 

A  operational  definition  of  Q  consistent  with  its  interpretation  as  a  triaxiality  pa¬ 
rameter  is 

Qm  =  — — at  0  =  0,  r  =  2J/oq  (2.15) 
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where  Om  is  the  hydrostatic  stress.  We  have  calculated  Q  based  on  the  hoop  stress  (2.9 )  and 
the  mean  stress  (2. 15)  for  the  full  range  of  T-stresses  and  finite  wdth  geometries.  The  dif¬ 
ference  between  Q  and  Qm  is  always  less  than  0.1. 

Our  numerical  studies  show  that  (2.14,2.15)  provide  a  more  accurate  representation  of 
the  full  range  of  near-tip  fields  so  that  a  fracture  methodology  based  on  (2.14,2.15)  has  a 
greater  range  of  validity.  Q  values  presented  in  this  paper  are  based  on  the  definition  in 
(2.9).  lb  simplify  subsequent  discussions,  the  term  small-scale  yielding  (SSY)  will  refer  to 
the  reference  field  (tJy )sSY;T*o  or  equivalently  (Oy )sSY;Q=0-  Fields  corresponding  to  IT^  0  or 
Q  0  will  be  explicitly  stated  as  such. 


2.6  Difference  Field  and  Higher-Order  Terms  of  the  Asymptotic  Series 


The  connection  between  the  difference  field  and  higher-order  terms  of  the  asymptotic  se¬ 
ries  pan  be  understood  in  the  context  of  the  MBL  formulation.  Here  the  stress  fields  obey 
the  functional  form 


(2.16) 


which  also  should  apply  to  finite— width  crack  geometries  as  long  as  the  characteristic  crack 
dimension  L  is  sufficiently  large  compared  taJloQ.  Now,  if  one  assumes  a  product  depen¬ 
dence  on  the  first  argument  in  (2.16)  and  works  with  deformation  plasticity  theory  with 
power  law  hardening  behavior,  then  one  obtains  a  series  in  r/(J  Ioq): 


<^0 


l/(n  +  l) 

ay(d,n)  +  second-order  terms  +  higher-order  terms 


(2.17) 


where  sq  is  a  reference  strain,  a  a  material  constant  (equal  to  unity  for  the  stress-strain 
response  defined  by  (2.10))  and  In  is  an  integration  constant.  By  definition,  the  asymptotic 
series  beyond  the  first  term  is  equivalent  to  the  difference  field  since  (see  Section  2.2  ) 


^ij  ~  difference  field  .  (2.18) 

Thus  the  HRR  field  and  (only)  the  second-order  term  provide  a  two-term  approximation 
to  the  solution  for  the  MBL  problem  and  this  point  appears  not  always  to  be  understood. 


The  higher  order  asymptotic  analysis  of  Li  and  Wang  [17]  and  Sharma  and  Aravas  [  18] 
has  been  extended  by  ^a,  Wang  and  Shih  [19].  They  have  obtained  a  five  term  expansion 
for  the  series  in  (2.17)  forn  =  3  and  a  four  term  series  for  n=10.  The  four  term  series  accu¬ 
rately  matches  the  radial  and  angxilar  variations  of  the  difference  field  given  in  Fig  3  and 
Fig.  5  by  O’Dowd  and  Shih  [13]  for  an  n  =10  material.  Indeed,  in  the  forward  sec¬ 
tor  1 0 1  <  jr/2,  the  collective  behavior  of  the  second,  third  and  fourth  order  terms  is  effectively 
equivalent  to  a  spatially  uniform  hydrostatic  stress  state  so  that  (2.17)  can  be  approxi¬ 
mated  by  the  simpler  form  in  (2.13). 


Finally,  it  may  be  noted  that  an  admissible  range  of  stress  states  for  an  elastic-perfectly 
plastic  material  can  be  written  in  the  form 

“  ^^y^Prandtl  »  1^1  ^  ^/4  ,  (2.19) 

where(ay)prandlt  designates  the  Prandtl  slip-line  solution  and  the  difference  field  corre¬ 
sponds  simply  to  a  iiniform  hydrostatic  stress  state  scaled  by  Q  [14,25]. 
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2.7  J-Q  Material  Toughness  Locus 


The  J-Q  theory  provides  the  quantitative  framework  to  characterize  a  material’s  fracture 
resistance  over  a  range  of  crack-tip  stress  triaxiality.  The  experimental  determination  of 
the  toughness  locus  and  its  utilization  in  engineering  applications  are  discussed  below. 

The  competition  between  fractiure  by  cleavage  and  ductile  tearing  controls  the  fracture 
resistance  of  ferritic  steels  in  the  ductile-to-brittle  transition  region.  Now  consider  test 
conditions  where  both  mechanisms  are  operative.  Fracture  by  (stress-controlled)  cleavage 
generally  requires  higher  crack-tip  constraint  while  ductile  tearing  develops  at  low 
constraint;  this  is  illustrated  by  the  two  distinct  segments  to  the  toughness  locii  shown  in 
Fig.  2.3a.  Since  measured  toughness  values  generally  exhibit  scatter,  both  the  lower  and 
higher  toughness  locii  are  indicated  which  define  bands  for  brittle  and  ductile  failure. 
Toughness  values  over  the  full  range  of  crack— tip  constraints  can  be  measured  by  using  the 
test  geometries  depicted  in  Fig.  2.3a.  As  an  example,  deeply-cracked  SE(B)  specimens  gen¬ 
erate  high  crack-tip  stress  triaxiality,  i.e.,  Q=0.  They  produce  driving  force  curves  which 
rise  steeply  and  therefore  intersect  the  toughness  locii  within  a  well-defined,  narrow  zone 
of  the  diagram.  In  contrast,  center-cracked  panels  and  single-edge  cracked  panels 
loaded  in  tension  are  low  constraint  crack  geometries.  'They  produce  driving  force  curves 
which  rise  with  more  shallow  slopes  and  thus  intersect  the  toughness  locii  over  a  broad  zone 
in  the  J-Q  diagram.  The  shallow  driving  force  curves  of  low  constraint  geometries  imply 
considerably  greater  scatter  in  cleavage  toughness  values  (Je),  a  phenomenon  commonly 
observed  in  testing  SE(B)  specimens  with  small  a/W  ratios,  for  example. 


(a)  Laboratory  Testing  (b)  Fracture  Assessment 

Figure  2.3  Application  of  the  J-Q  methodology  in  fracture  assessments,  (a)  Laboratory  test¬ 

ing  of  specimens  with  varying  constraint  to  measure  the  material’s  fracture  resis¬ 
tance.  Circles  indicate  anticipated  scatter  which  define  upper-lower  bounds,  (b) 
Evaluation  of  structural  flaws  using  measured  toughness  locus  and  predicted  J-Q 
response  for  two  structural  configurations.  Cleavage  fracture  is  predicted  for 
Structure  A;  ductile  tearing  is  predicted  for  Structure  B. 

Utilization  of  the  toughness  locus  in  fracture  assessments  is  illustrated  in  Fig.  2.3b. 
Suppose  that  the  material’s  fracture  resistance  under  service  conditions  is  characterized 
by  the  indicated  cleavage-ductile  failiire  band.  The  driving  force  curve  for  a  structure  of 
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high  crack-tip  constraint,  structure  A,  rises  rapidly  in  the  J-Q  space  so  that  cleavage 
fracture  occurs  when  the  driving  force  curve  intersects  the  fe’lure  locus.  In  contrast  a  low 
constreiint  geometry  such  as  structvu-e  B,  induces  a  gradually  rising  driving  force  curve  so 
that  ductile  tearing  is  the  likely  event  at  overload. 
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3.  MICROMECHANICAL  CONSTRAINT  CORRECTIONS 

Dodds  and  Anderson  [31,32]  show  that  by  quantifying  the  effects  of  finite  size  on  the  rela¬ 
tionship  between  microscale  crack  driving  force  (e.g.  near-tip  stresses  and  strains)  and 
macro-scale  crack  driving  force  (e.g.  J,  CTOD),  the  apparent  size  effect  on  fracture  tough¬ 
ness  can  be  predicted  rigorously  without  resort  to  empiricad  arguments.  These  size  effects 
become  steadily  more  pronounced  as  load  increases  due  to  the  deviation  of  crack-tip  region 
deformations  from  the  small  scale  )delding  conditions  essential  for  single  parameter  frac¬ 
ture  mechanics  (SPFM)  to  apply.  VWien  SPFM  becomes  invalid,  a  micromechanics  failure 
criteria  is  required  to  establish  the  near-tip  conditions  at  fracture.  Finite  element  analysis, 
or  alternatively  the  near-tip  stresses  described  by  the  J-Q  theory,  provides  a  means  to 
quantify  the  geometry  dependent  relations  between  these  micromechanical  failure  condi¬ 
tions  and  macro-scale  crack  driving  force.  This  permits  (in  principle)  prediction  of  fracture 
in  any  body  from  toughness  values  measured  using  standard  specimens. 

For  steels  operating  at  temperatures  where  cleavage  occurs  after  significzint  plastic  de¬ 
formation  but  before  the  initiation  of  ductile  growth  (lower  to  mid-transition),  attainment 
of  a  critical  stress  over  a  microstructursdly  relevant  volume  is  an  appropriate  microme¬ 
chanical  failure  criteria  [21-39].  A  niimber  of  important  engineering  structures  can  fail  by 
this  mechanism,  including  high  strength  rails,  offshore  oil  platforms,  ships,  storage  tanks, 
and  nuclear  pressure  vessels  after  years  of  neutron  irradiation  embrittlement.  Techniques 
for  predicting  the  apparent  size  effects  on  cleavage  fracture  toughness  developed  by  Dodds 
and  Anderson  are  described  in  the  following  sections. 

3.1  Transgranular  Cleavage  Mechanism 

A  number  of  micromechanical  models  for  transgranular  cleavage  fracture  have  been  pro¬ 
posed,  most  derive  from  weakest-link  statistics.  The  weakest-link  models  assume  the  Isu'g- 
est  or  most  favorably  oriented  fractvure-triggering  particle  controls  the  cleavage  failure. 
The  actual  trigger  event  involves  a  local  Griffith  instability  of  a  microcrack  which  forms  at 
a  microstructural  feature  such  as  a  carbide  or  inclusion;  satisfaction  of  of  the  Griffith  ener¬ 
gy  balance  occurs  when  the  critical  stress  is  reached  in  the  vicinity  of  the  microcrack.  The 
size  and  location  of  the  triggering  microstructvural  feature(s)  dictate  the  fracture  toughness 
and  produces  the  scatter  routinely  observed  in  results  of  cleavage  fracture  tests. 

The  Griffith  instability  criterion  implies  fracture  at  a  critical  normal  stress  near  the 
crack  tip;  the  statistical  sampling  aspect  of  the  mechanism  (i.e.,  the  probability  of  finding 
a  triggering  microfeature  near  the  crack  tip)  suggests  a  dominant  role  for  the  volume  of  ma¬ 
terial  within  a  process-zone  over  which  the  opening  mode  stress  exceeds  a  threshold  value 
sxifficient  to  initiate  cleavage.  The  probability  of  cleavage  fracture  in  a  cracked  specimen 
may  then  be  expressed  in  the  following  general  form: 

F  *  FTV(ai)]  (3.1) 

where  F  is  the  failure  probability,  oi  is  the  maximum  principal  stress  at  a  point  and  V(ai) 
is  the  cumulative  voltime  sampled  over  which  the  principal  stress  is  equal  to  or  greater  than 
ai.  This  form  of  F  applies  to  any  fracture  process  controlled  by  maximum  principal  stress, 
not  just  weakest-link  failure  which  is  now  being  questioned  [40,41].  In  particular,  the  F 
criterion  of  (3.1)  does  not  require  material-specific  assumptions  for  the  distribution  and 
strength  of  cleavage  triggering  particles. 

Unlike  other  micromechanics  models,  the  present  methodology  does  not  attempt  to  pre¬ 
dict  absolute  values  of  Jc  from  metallurgical  parameters  that  describe  the  distribution  and 
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Strength  of  cleavage  triggering  paulicles.  Rather,  the  micromechanical  model  predicts  the 
variation  of  fracture  toughness  with  constraint  changes  for  a  given  materiaL/temperature 
by  scaling  to  a  reference  condition.  The  crack-tip  stress  fields  in  a  test  specimen  are 
compared  to  the  limiting  solution  of  SSY.  A  J-like  parameter,  denoted  Jq.  is  obtained  from 
this  comparison  to  the  reference  solution.  Jo  'is  the  J  to  which  the  SSY  model  (infinite  body) 
must  be  loaded  to  achieve  the  same  stressed  volume,  and  thereby  the  same  likehhood  of 
cleavage  fracture,  as  in  a  finite  body. 

A  critical  value  of  Jo  represents  the  fracture  toughness  of  an  infinitely  large  specimen; 
the  ratio  of  applied  J/Jo>  1  implies  that  the  specimen  has  experienced  a  constraint  loss 
that  causes  the  commonly  observed  increase  in  measured  fractured  toughness. 

3.2  Development  of  Constraint  Corrections 


By  employing  the  family  of  near-tip  states  in  the  form  of  (2.4),  the  maximum  principal 
stress  also  has  the  form 


‘'J  ■ 


For  any  given  value  of  Q  and  6,  OiIoq  decreases  monotonically  once  r  extends  beyond  the 
finitely  deformed  region  of  r<J Ioq,.  Rearrangement  of  the  above  expression  furnishes  a 
relation  for  the  distance  r  as  a  fimction  of  Q  and  ai/ao  as 

r  =  ^gi(0;£ri/ao,Q)  .  (3.3) 

Consider  a  particular  level  of  the  principal  stressai/ao.  TheareaA  over  which  the  principal 
stress  is  greater  than  cri/ao  is  given  by 


A  =  Q),  ^  “  i  I  oJoq,  Q)d0 

^0  J 


The  area  enclosed  by  the  contour  of  level  oilao  depends  on  J  as  well  as  the  triaxiality  of 
the  near-tip  fields  identified  with  Q.  Tb  fix  ideas,  let  Aq  and  Jo  designate  the  area  and  J 
eissociated  with  the  Q=0  field,  and  let  Afb  aiid  Jfb  designate  the  area  and  J  Eissociated  with 
a  crack  in  a  finite  body  with  Qf^O.  Then  we  have 

tP  f 

^0  =  ^0  =  ^1  Q  =  O)d0  (3.5) 


^  r 

•^FB  “  ^FB  ~  ^  \  •  (3.6) 

Upon  initial  loading  of  the  finite  body,  Q=0  so  that  hj^=ho;  compare  expressions  (3.5b)  and 
(3.6b).  As  the  load  increases,  plasticity  spreads  over  the  body,  Q  becomes  non-zero,  and  hpB 
begins  to  deviate  from  ho. 

For  a  given  material  and  temperature,  the  present  micromechanics  model  requires  the 
attainment  of  equivalent  stressed  volumes  (ApB  x  thickness  B)  for  cleavage  fracture  in  dif¬ 
ferent  specimens.  The  ratio  of  applied  J-values  in  a  finite  body  and  the  reference  Q=0  stress 
state  that  generate  equivalent  stressed  volumes  is  found  by  equating  areas  in  (3.5)  and 
(3.6)  to  yield 
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«^FB  _  /  hQ{ai/a^  ^2 

Jo  V 

The  resvJts  thus  far  are  generally  applicable  and  do  not  rely  upon  any  particular  form  of 
the  J-Q  fields. 

The  J  ratios  are  evaluated  using  (3.7)  at  each  loading  level  and  for  a  range  of  principal 
stress  values.  The  ratio  quantifies  the  size  and  geometry  dependence  of  cleavage  fracture 
toughness.  Consider,  for  example,  a  test  specimen  that  fails  at  t7c=200  kPa  •  m.  Suppose  the 
computed  ratio  Jfb/‘^o=2  at  fracture  (Jfb=^c)  ^  test  specimen;  then  a  very  large  speci¬ 
men  made  firom  the  same  material  and  tested  at  the  same  temperature  is  predicted  to  fail 
at  Jc=100  kPa  •  m.  Similarly,  the  fracture  toughness  ratios  for  test  specimens  with  the  same 
absolute  size  but  varying  crack— depths  to  specimen-widths,  a/W,  may  be  quantified.  The 
model  predicts  a  sharp  increase  in  fracture  toughness  with  decreasing  a /IV  ratio. 


Self-Similar  Principal  Stress  Contours 


The  character  of  the  near-tip  fields  has  been  investigated  by  O’Dowd  and  Shih  [13,14]  and 
Xia,  Wang  and  Shih  [19].  From  their  results,  e.g.,  (2.14),  we  write 


(3.8a) 

(3.8b) 


where  the  form  in  (3.8b)  describes  the  fields  in  the  forward  sector,  \0\<7tl2  and  r<5J  I  oq. 
Detailed  computational  studies  have  shown  that  principal  stresses  of  s\ifriciently  high  lev¬ 
el,  say  oi/oq  >  2.0,  are  found  only  in  the  forward  sector.  The  form  in  (3.8b)  is  applicable  in 
the  preceding  micromechanics  analysis  if  we  confine  attention  to  ai  Ioq  >  2.0.  Rearranging 
(3.8b)  yields 


Zl 

^0 


(3.9) 


The  form  in  (3.9)  constitutes  a  self-similar  field  for  cti/oq-  Q.  Moreover,  the  behavior  of 
oi/oq-Q  obeys  the  form  governing  the  Q=0  reference  solution,  which  has  been  determined 
by  small-scale  yielding  analysis. 


lb  understand  the  implications  of  (3.9),  we  focxis  on  a  particular  value  of  ailao,  say 
ai/a(y=3.  Consider  the  behavior  of  oi/ao  -  Q  as  the  deformation  level,  measiired  by  J,  in¬ 
creases.  Suppose  for  the  moment  that  Q  remains  constant;  then  the  contour  for  a  fixed  level 
of  ai /ao- Q,  presented  in  the  normalized  distancesX/(J /oq)  and  Y/{J /oq),  remains  unal¬ 
tered  in  size  and  shape  with  increasing  deformation  level.  As  an  example,  the  outermost 
contour  in  Fig.  3.1  corresponding  to  the  reference  solution,  fixed  Q  =  0,  maintains  its  size 
and  shape  as  the  deformation  level  is  increased. 

Now  let  Q  evolve  with  increasing  deformation  as  happens  in  a  finite  size  body.  Q  de¬ 
creases  gradually  corresponding  to  a  loss  of  stress  triaxiality  as  the  deformation  level  in¬ 
creases;  oi/ao-Q  must  increase  since  a\/ao  is  fixed.  Therefore  the  evolution  (size  reduc¬ 
tion)  of  a  contour,  associated  with  a  fixed  value  of  oi  /oq,  vmder  increasing  plastic  yielding 
can  be  described  by  contours,  associated  with  increasingly  higher  levels  of  oiIoq- Q,  gov¬ 
erned  by  (3.9).  The  sequence  of  diminishing  contours  associated  with  increasing  levels  of 
/OQ-  Q,  corresponding  to  a  fixed  level  of  ff  i  / ao=3,  is  depicted  in  Fig.  3. 1  for  a  shallow  notch 
SE(B)  specimen. 
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Figure  3.1  Comparison  of  a  maximum  principal  stress  contour  for  SSYq=o  with  those  for 
an  o/ W=0.15,  n=10  SE(B).  SE(B)  contours  decrease  in  size  with  increasing  de¬ 
formation  (i.e.  with  decreasing  o(7o/J). 


2.0  2.5  3.0  3.5  4.0 


Figure  3.2  Areas  within  principal  stress  contours  for  an  a  /  W=0. 15,  n- 10  SE(B).  Values  are 

normalized  by  area  within  contour  for  SSYqso  at  same  J-value. 
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Finally,  the  asymptotic  studies  show /o  is  nearly  separable  in  r  and  Q  when  r/(J/ao)  is 
sufficiently  small  (HRR  field).  Consequently,  the  shape  of  a  contour  is  maintained  as  its  size 
diminishes.  The  self-similarity  of  oiIoq-Q  prevails  only  to  the  extent  that  the  J-Q  form 
in  (2.14)  remains  applicable. 

Inverted  Relationships 


A  specialized  form  of  (3.3)  is  developed  by  noting  the  dependence  of  r  on  ai  /ao  and  Q  involve 
them  in  the  combination  oiIoq  -  Q,  i.e., 

r  =  4-^gT(.6-,aJa^  -  Q)  .  (3.10) 

The  results  in  (3.4)  through  (3.7)  are  simplified  by  using  the  form  in  (3.10). 

The  key  question  to  resolve  with  this  approach  concerns  the  sensitivity  of  the  Ayr/Aq 
and  ratios  to  OiIoq.  Let  Jj>b  denote  the  value  associated  with  (ai  /ctq)*.  Then  to  first- 

order. 


When  Q  *=  0,  the  ratio  JpB  Uo  is  insensitive  to  txi / ao  since  the  quantity  in  []  is  scaled  by  Q. 
When  Q  is  large  (negative),  the  JpB  Uo  exhibits  a  small  sensitivity  to  ai  / ao.  The  weak  de¬ 
pendence  of  Jfb/«/o  on  oiloQ  has  been  confirmed  by  analyzing  the  evolution  of  near-tip 
fields  in  common  fracture  specimens.  Figures  3. 1-3.3  provide  typical  results  obtained 
through  finite  element  modeling.  The  specimen  is  a  single-edge  notched  bend  bar  contain¬ 
ing  a  shallow  notch,  a/W=  0. 15,  with  a  strain  hardening  exponent  of  a=  10.  Figure  3.2  shows 
the  area  enclosed  by  principal  stress  contours  (Afb)-  The  SE(B)  areas  are  normalized  by  the 
area,  Ao,  defined  by  the  same  contour  of  the  reference  solution  (Q=0)  when  loaded  to  the 
same  J  as  the  SE(B);  Jo=Jyb-  The  area  ratios  remain  relatively  insensitive  to  oi/oq  until 
the  deformations  become  excessive.  «7o  is  calculated  for  each  line  of  this  figure  using  (3.7). 
Jy^Jo  ratios  are  independent  of  the  principal  stress  selected  for  computation  over  a  wide 
range  as  shown  in  Fig.  3.3.  In  practice,  the  computation  of  t/pB/t/o  ratios  is  terminated  when 
the  values  differ  by  more  than  10%  at  the  smallest  and  largest  principal  stress  values  as 
indicated  on  the  figure.  A  larger  deviation  indicates  breakdown  in  the  similarity  of  the 
SE(B)  and  SSY  stress  fields  and  thus  an  unacceptably  large  dependence  on  the  critical  frac¬ 
ture  stress. 


To  simplify  applications  of  this  methodology,  the  SSY  areas  within  principal  stress  con¬ 
tours  are  expressed  as  a  function  of  principal  stress  (ai)  and  strain  hardening  coefficient 
in)  as: 


-  -  =  10^/250000 

■'o 

where  the  curve  fitting  function  is  given  by 


(3.12) 


with  fitting  coefficients  Hi  given  in  Table  3.1  for  a  range  of  hardening  exponents. 


(3.13) 
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Figure  3.3  Influence  of  specified  critical  stress  on  the  micromechanics  prediction  of  frac¬ 
ture  toughness  variation  with  constraint  for  an  o/W=0.15,  n=10  SE(B). 
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Figure  3.4  Micromechanics  prediction  of  alW  effects  on  cleavage  fracture  toughness  for  an 
n=5,  SE(B)  fracture  specimen. 
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Table  3.1:  Fit  coefficients  for  (3.13) 


n 

Ho 

Hi 

■  1 

H2 

Hs 

H4 

Mini¬ 

mum 

oi/ao 

Maxi¬ 

mum 

Oi/Oo 

4.0 

6.4306 

-2.4711 

0.5037 

-0.07975 

0.00552 

2.0 

4.0 

5.0 

6.2579 

-2.1653 

0.3749 

-0.06603 

0.00505 

2.0 

4.0 

10.0 

7.6641 

-4.3138 

1.7368 

-0.43685 

0.03560 

2.0 

3.6 

20.0 

-3.2613 

14.4338 

-10.2659 

3.01033 

-0.34420 

2.0 

3.2 

These  expressions  are  obtained  by  curve  fitting  the  results  of  small-displacement  gradient, 
finite  element  analyses  conducted  on  the  MBL  problem  with  T/ao=0.  The  material  stress- 
strain  curve  employed  in  the  analyses  follows  the  conventional  Ramberg-Osgood  model 
given  by 


which  exhibits  a  slightly  different  behavior  than  the  elastic  power-law  model  defined  in 

(2.10). 

3.3  Application  of  the  Constraint  Corrections  in  Fracture  Testing 

The  computational  procedures  outlined  above  have  been  applied  to  generate  t/pB/^o  ratios 
for  a  variety  of  test  specimens  and  material  properties  [30,31].  Figure  3.4  provides  the  re¬ 
sults  of  such  computations  for  SE(B)  specimens  having  a  range  ofa/W  ratios  modelled  with 
an  n-5  strain  hardening  material.  Values  of  JpB  and«7o  are  plotted  on  sepai  Ate  axes  to  facil¬ 
itate  removal  of  the  size  effect  in  experimental  data.  Points  on  the  curves  describe  ( JpBi  Jo) 
pairs  that  produce  equal  stressed  volumes  of  material  in  the  finite-size  test  specimen  and 
in  the  SSY  model.  Upon  initial  loading,  crack-tip  pleisticity  is  well  contained  within  a  sur¬ 
rounding  elastic  field  and  identical  values  for  JpB  and  Jq  correspond  to  the  same  stressed 
voluiiie  of  material  at  the  crack  tip.  This  1:1  line  is  shown  on  the  figure  for  reference.  At 
higher  loads  and  as  constraint  relaxes  under  extensive  plastic  flow,  the  finite-size  test  spec¬ 
imen  requires  more  applied- J  («/fb  >«^o)fo  achieve  the  SEune  conditions  for  cleavage  (same 
stressed  volume)  as  in  SSY. 

Information  of  this  type  is  useful  for  both  analysis  of  fracture  test  data  and  for  assessing 
the  defect  integrity  of  structures.  Path  A-B-C  on  Fig.  3.4  illustrates  the  procedure  to  re¬ 
move  geometry  dependence  from  experimental  cleavage  fracture  toughness  data  (Jc value 
at  A)  by  determining  the  geometry  independent  cleavage  fracture  toughness  (Jo  value  at 
C)  corresponding  to  a  measured  Jc  value.  Alternatively,  Fig.  3.4  permits  determination  of 
the  apparent  fracture  toughness  for  an  SE(B)  with  any  a/W  ratio  firom  a  known,  Jq  value 
(path  C-D-E  for  example). 

Figure  3.5a  shows  the  Jc  values  measured  by  Sumpter  and  Forbes  [10]  for  a  BS4360 
43A  steel  in  =>5)  using  SE(B)  specimens  tested  over  a  wide  range  of  a/W  ratios.  The  data 
readily  demonstrate  the  dramatic  increase  in  cleavage  fracture  toughness  with  decreasing 
a/W  ratio.  To  remove  the  constraint  effect  on  toughness,  each  experimental  data  point  is 
processed  using  a  path  similar  to  A-B-C  in  Fig.  3.4  to  obtain  the  corresponding  Jq  value. 
Figure  3.5b  shows  these  “constraint  corrected”  toughness  values.  The  toughness  variation 
with  a/W  ratio  is  effectively  removed  with  this  technique.  The  small  remaining  scatter  in 
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the  experimental  data  for  different  aAV  ratios  may  be  attributed  to  true  metallurgical  vari¬ 
ations  in  the  material  and  the  unavoidable  procedural  variations  in  testing  a  large  number 
of  specimens.  Application  of  this  technique  to  other  materials  including  A36,  A515  and 
A533B  [42]  have  been  equally  successful  in  removing  the  geometry  dependence  of  Jc  val¬ 
ues. 

Jc  [in  •  Ibs/in^]  Jq  Lin  •  lbs /in}] 


Figure  3.5  (a)  Effect  of  initial  cracV  depth  on  cleavage  fracture  toughness  in  a  mild  steel 

(Sumpter  and  Forbes  [10]);  (b)  Jo  values  (specimen  size  independent  fracture 
toughness)  calculated  from  experimental  Jg  data  using  the  micromechanics 
constraint  coi  rection. 

3.4  Engineering  Use  of  J-Q  Fields  in  the  Micromechanics  Model 


When  available,  a  J-Q  description  of  the  crack-tip  stresses  for  a  test  specimen  or  structural 
component  may  be  readily  employed  to  generate  the  constraint  corrections  for  fracture 
toughness  of  the  type  shown  in  Fig.  3.4.  Here  we  outline  a  proced\ire  that  is  computation¬ 
ally  simpler  than  the  stressed  volume  approach  defined  by  (3.2)-(3.7)  but  which  yields  es¬ 
sentially  the  same  result. 


Figure  3.6  shows  the  variation  of  opening  mode  stress  on  the  crack  plane  with  distance 
from  the  crack  tip  for  several  deformation  levels  for  an  SE(B)  specimen  with  a/W=0.15, 
n=10.  The  SE(B)  stresses  are  normalized  by  the  stress  in  the  SSY  model  at  the  same  rela¬ 
tive  distance  ahead  of  the  crack  tip  when  the  SSY  model  is  loaded  to  the  same  J  as  the 
SE(B).  Distances  are  normalized  by  the  similarity  length-scale  r/{J/{aafff)).  The  indepen¬ 
dence  of  these  normalized  stresses  with  distance  from  the  crack  tip  indicates  again  the  simi¬ 
larity  of  the  SSY  and  SE(B)  stress  distributions.  Jq  is  calculated  at  a  number  of  points  along 
each  line  on  this  graph  as  the  J  value  required  in  the  SSY  model  to  achieve  the  same  open¬ 
ing  mode  stress  as  in  the  finite  body.  The  following  equation  is  solved  iteratively  for  Jq  using 
a  simple  nonlinear  root  solver: 


at  0  =  0 

where  a  closed-form  fit  to  the  crack-plane  stresses  in  SSY  is  given  by 


(3.15) 


where 


(3.16) 
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stress  at  distance  r  from  the  tip  at  loading  Jfb  in  finite  body; 


r  normalized  distance  from  tip:  r  =  -  - ; 

Gj  SSY  fit  coefficients  summarized  in  Table  3.2. 

The  finite-body  stresses  needed  in  (3.15)  are  given  by 


•<»)pB  _  (o'eehsY  , 

^0 - 


(3.17) 


where  the  SSY  term  in  (3.17)  is  given  by  the  expression  on  the  right  side  of  (3.16)  with  Jq 
replaced  by  JpB-  The  potential  for  a  small  radial  dependence  of  Q  under  large-scale  yielding 
is  included  in  (3.17).  The  JpB/efo  ratio  is  computed  over  a  range  of  distances  ahead  of  the 
tip  at  each  load  level,  typically  1-2  s  r/iJ/ao)  s  4-5  with  the  specific  values  dependent  upon 
the  degree  of  strain  hardening.  The  objective  is  to  sample  the  stress  field  at  locations  out¬ 
side  the  finitely  deformed  zone,  r>2d,  but  within  the  process  zone  applicable  for  cleavage 
fractinre.  Herrens  and  Read  [43];  Miglin,  et  al.  [44]  determined  fractographically  the  limit 
of  the  cleavage  process  zone  as  r>s8<5. 


Figure  3.7  shows  the  Jfb/Jq  ratios  computed  using  this  approach  for  each  loading  level 
indicated  in  Fig.  3.6.  The  similarity  between  the  SSY  and  finite-body  stress  distributions 
makes  the  specific  r/d  value  used  in  the  calculations  unimportant  over  a  wide  range  of  de¬ 
formation.  In  practice,  Jq  calculated  by  (3.15)  is  considered  valid  when  the  values  calcu¬ 
lated  at  r=3d  and  at  r=8d  differ  by  less  than  10%.  A  larger  deviation  signals  too  great  a 
dependence  of  Jq  on  the  critical  distance  selected  and,  consequently,  a  breakdown  of  the 
method.  Figure  3.8  compares  the  constraint  corrections  for  fracture  toughness  computed 
using  the  simpler  approach  with  crack-plane  stresses  given  by  a  J-Q  analysis  and  the  more 
complex  approach  requiring  computation  of  stressed  volumes  within  principal  stress  con- 
toms.  Differences  in  the  constraint  corrections  are  insignificant  for  engineering  applica¬ 
tions.. 


n 

4.0 

5.0 

7.0 

Table  3J2:  Fit  coefficients  for  SSYqbO  iu  (3.16) 
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-0.1169 
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18.0 
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Figure  3.6  Opening  stress  on  the  crack  plane  normalized  by  SSYqeO  at  same  (/-applied. 
a/VY=0.15,  n=10  SE(B). 


r/d 


Figure  3.7  Variation  of  fracture  toughness  predicted  by  the  micromechanics  model  using 
crack-plane  stress  shown  in  Fig.  3.6.  a/lV=0.15,  n=10  SE(B). 
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Figure  3.8 
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Comparison  of  Jo  values  determined  using  the  crack-plane  stress  from 
J-Q  and  stressed  volume  techniques  for  an  a/ W=0.15  n=10  SE(B). 
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4.  SURFACE  CRACKS  UNDER  BIAXIAL  LOADING 


Bass,  et.  al  [45]  recently  outlined  current  deficiencies  in  the  understanding  of  constraint 
effects  on  the  crack-initiation  toughness  of  shallow  surface  cracks  subjected  to  uniaxial 
and  biaxial  far-field  tension  loadings.  In  nuclear  applications,  the  internal  pressure  alone 
generates  a  1:0.5  biaxial  tension  loading  while  the  pressurized  thermal  shock  (PTS)  event 
may  generate  the  more  severe  case  of  1:1  biaxial  tension  in  addition  to  a  locally  severe  bend¬ 
ing  field.  The  very  few  testing  programs  (see  Bass,  et.  al  [45]  for  a  review)  conducted  on 
biaxiedly  loaded  surface  cracks  report  a  25—40%  reduction  in  toughness  values  (XJ;)  relative 
to  the  values  obtained  from  SE(B)  and  C(T)  specimens  containing  cracks  of  similar  relative 
depth.  These  results  imply  a  significantly  increased  crack-tip  constraint  under  biaxial 
loading  relative  to  \iniaxial  loading.  Moreover,  the  biaxial  test  resiilts  appear  to  negate  the 
now  well  established  increases  in  fracture  toughness  for  shallow  notch,  SE(B)  specimens 
relative  to  deep  notch  toughness. 

Researchers  currently  frame  discussions  of  the  biaxial  vs.  uniaxial  loading  influence  on 
constraint  in  terms  of  in-plane  and  out-of-plane  effects.  Shallow  crack  SE(B)  specimens, 
for  example,  exhibit  a  strong  in-plane  effect  on  constraint;  the  small  crack  depth  relaxes 
crack-tip  stresses  when  plastic  zones  sense  the  nearby  fi’ee  siirfaces  behind  the  crack.  Out- 
of-plane  effects  refer  to  tensile  stresses  acting  parallel  to  the  crack  front.  While  these 
stresses  exist  and  vaiy  along  the  front  of  uixiaxially  loaded  specimens,  test  programs  dem¬ 
onstrate  the  much  smaller  influence  of  thickness  (E),  which  governs  the  out-of-plane 
stress,  relative  to  the  crack-depth  effect  once  B  exceeds  a  significant  fraction  of  the  speci¬ 
men  width  (W),  usually  B  >  W/2.  The  biaxial  test  results  suggest  that  mechanically  applied, 
remote  out-of-plane  stress  restores  crack-tip  triaxiality  lost  to  the  shallow-crack  in-plane 
effect.  Strength-of-materials  type  models  have  been  proposed  to  examine  the  interaction 
of  in-plane  and  out-of-plane  stresses.  Such  methods  are  severely  limited  since  they  rely 
on  superposition  of  stresses  which  does  not  apply  under  elastic-plastic  conditions  at  the 
crack  tip. 

The  scarcity  of  testing  programs  that  address  biaxial  loading  effects  on  fracture  tough¬ 
ness  and  their  significant  complexity  (large  plate  specimens,  scale-model  pressure  vessels, 
thermo-mechanical  loading,  etc.)  leave  open  an  experimental  resolution  of  this  issue.  How¬ 
ever,  the  J-Q  and  micromechanics  concepts  described  in  Sections  2  and  3  provide  the  ana¬ 
lytical  framework  to  clarify  the  in-plane  and  out-of-plane  effects  on  crack-tip  stress  fields 
(using  J-Q)  and  on  cleavage  fracture  toughness  (using  micromechanics). 

4.1  Part-Through  Surface  Crack  Model 

Figure  4.1  shows  a  flat  plate  containing  a  part-through  surface  crack  considered  in  an  ini¬ 
tial  analytical  investigation.  The  semi-elliptical  surface  crack  has  geometric  parameters 
a/f=0.25,  2c/as6.0.  <t>=Q  and  ^=90°  correspond  to  lines  along  the  free  surface  and  directly 
ahead  of  the  front  at  the  point  of  maximum  depth.  Radial  distances  ahead  of  (and  normal 
to)  the  crack  front  are  denoted  r.  The  plate  is  loaded  by  a  remote  uniaxial  tension,  oZ  in 
one  case,  and  by  a  remote  biaxial  tension,  aZ  =  in  the  second  case.  The  material  re¬ 
sponse  is  modeled  with  small-strain,  deformation  plasticity;  the  uniaxial  flow  properties 
obey  a  Ramberg-Osgood  relationship  with  hardening  exponent  n=10  and  a=l. 

S}naimetry  conditions  enable  consideration  of  only  one-quarter  of  the  full  specimen  in 
the  f^te-element  model  as  shown  if  Fig.  4.1.  The  element  mesh  contains  1980  20-node, 
isoparametric  elements  and  9800  nodes.  The  level  of  mesh  refinement  in  r  at  each  point 
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Figure  4.1  Finite-element  model  for  investigation  of  constraint  in  surface  cracked  plate 
subjected  to  uniaxial  and  biaxial  remote  tension  loadings. 
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Figure  4.2  Global  load-displacement  (Crack  Mouth  Opening  Displacement,  CMOD)  re¬ 
sponse  for  surface  cracked  plate. 
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Figure  4.3  J-integral  values  at  point  of  maximuin  crack  depth  for  uniaxial  and  biaxial  ten¬ 

sion  loading  of  siuface  cracked  plate  model. 


Figure  4.4  (/-integral  along  crack  front  at  SSY  and  LSY  load  levels  normalized  by  corre¬ 
sponding  J  at  maximum  crack  depth  for  each  loading  level. 
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along  the  crack  front  matches  the  refinement  employed  in  previous  models  [  14,31]  of  two- 
dimensional  specimens.  The  innermost  ring  of  elements  incident  on  the  crack  front  con¬ 
tains  degenerate  20-node  elements  with  edge  nodes  retained  in  the  mid-point  location.  Ini- 
tieiUy  coincident  nodes  along  the  crack  front  are  unconstrained  to  permit  blunting  deforma¬ 
tions.  Uniform  reduced  integration  (2x2x2)  in  all  elements  performed  satisfactorily  in  these 
models. 


The  intensity  of  local  deformation  at  each  point  s  along  the  front  is  given  by  [46] 


du. 


dr 


(4.1) 


where,  W  denotes  the  strain-energy  density,  7^  is  a  vanishingly  small  contour  in  the  princi¬ 
pal  normal  plane  at  s,  n  is  a  unit  normal  vector  to  7^,  Oy  and  uj  are  Cartesian  components 
of  stress  in  the  crack  front  coordinate  system.  Numerical  evaluation  of  (4. 1)  is  accomplished 
with  a  domain  integral  method  [46,47]. 

Figure  4.2  shows  the  overall  load-displacement  response  in  terms  of  Crack  Mouth 
Opening  Displacement  (CMOD).  Under  SSY,  CMOD  remains  unaffected  by  the  biaxial 
loading.  With  the  onset  of  gross  plasticity,  however,  the  biaxial  loading  provides  a  signifi¬ 
cant  stiffening  effect;  at  (7“/c7o=l.l  the  biaxial  CMOD  is  only  55%  of  the  uniaxial  value. 


A  similar  effect  of  the  biaxial  loading  on  theJ-values  can  be  seen  in  Fig.  4.3.  At  the  point 
of  maximum  crack  depth  (^=90°),  the  uniaxial  J-value  is  twice  the  biaxial  value  when  both 
models  are  loaded  to  a*/(7o=l.l.  The  comparison  of  biaxial  and  uniaxial  distributions  for 
Jlocal  along  the  crack  front  is  shown  in  Fig.  4.4.  The  distributions  are  identical  under  SSY 
but  reveal  considerable  differences  under  large  scale  yielding  in  the  region  of  sharpest  front 
curvature  (4>  <  30°).  The  biaxial  loading  depresses  the  level  of  J  relative  to  the  uniaxial  load¬ 
ing  in  this  region  of  the  crack  front. 


4.2  Crack-Front  Stress  Triaxiality 


Figure  4.5  shows  the  behavior  of  near-tip  stress  triaxiality  along  radial  lines  normal  to  the 
crack  front  at  ^=17°  and  90°  for  the  uniaxial  and  biaxial  loadings.  Q-values  are  defined 
from  opening-mode  stresses  on  the  crack-plane  using  (2.9).  At  <p=  90°,  the  variation  of  Q 
with  r  remains  negligible  up  to  the  maximum  applied  load  of  a“/ao=l.l.  At  loadings 
a"/aoS0.4,  the  model  lacks  sufficient  refinement  to  resolve  stresses  over  the  region 
2  ^  ^  5.  The  crack-tip  constraint  steadily  decreases  with  increased  global 

loading  and  plastic  deformation  (Q  becomes  more  negative).  The  biaxial  loading  exerts  only 
a  minor  influence  toward  reducing  the  constraint  loss  under  large-scale  yielding  at  <p=  90°. 

The  ^=17°  location  on  the  crack  front  has  high  cxirvature  and  lies  a  smsdl  distance  from 
the  traction-free  face  of  the  plate  (see  Fig.  4.1).  Under  increased  uniaxial  loading,  Q  steadi¬ 
ly  decreases  indicating  a  gradual  loss  of  stress  triaxiality.  Q  develops  a  weak  dependence 
on  radial  distance  as  the  maximum  applied  loading  is  approached;  Q  varies  by  ±  6.7%  from 
the  mean  value  over  2  ^  ^  5  at  a*/ao=l.l.  The  biaxial  loading  maintains 

stress  triaxiality  at  significantly  higher  levels  as  plastic  flow  progresses  from  well-con¬ 
tained  through  large-scale  yielding.  The  radial  dependence  increases  to  ±17%  over 
2  ^  s  5  at  a*/ao=l.l.  The  radial  dependence  of  Q  when  large-scale  yielding 

prevails  is  induced,  in  large  part,  by  the  nearby  free  surface. 

Alternatively,  Q  may  be  computed  from  values  of  the  mean  stress  ahead  of  the  crack  tip 
by  appealing  to  the  defining  equation  (2.15).  Figure  4.6  illustrates  the  potential  differences 
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(a)  Uniaxial  Tension 


loctdf^O^  local/ 

(b)  Biaxial  Tension 


Figure  4.5  Evolution  of  stress  triaxiality  in  surface  cracked  plate  ahead  of  crack  tip  at  ^=90' 

and  17®. 
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Figure  4.6  Comparison  of  Q-values  computed  using  opening-mode  stress  (2.9)  and  the 
mean  stress  (2.14):  (a)  SE(B)  specimen  with  o/W=0.15,  n=10;  (b)  Surface- 
cracked  plate  shown  in  Fig.  4.1. 
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in  the  two  methods  to  evaluate  Q  (2.9,  2.15).  Q-values  computed  from  the  opening  mode 
stress  and  the  mean  stress  are  showm  for  a  plane-strain  model  of  a  shallow  notch  SE(B) 
specimen  and  for  the  present  surface  crack  model.  The  SE(B)  specimen  exhibits  no  differ¬ 
ence  for  the  two  computational  procedures  while  the  small  differences  at  high  load  levels 
for  the  surface  crack  are  considered  insignificant  for  engineering  applications.  The  Q-val- 
ues  for  the  surface  crack  shown  in  Fig.  4.6  represent  the  condition  (uniaxial  vs.  biaxial  load¬ 
ing,  crack-front  location)  having  the  largest  discrepancy  for  the  two  computationcil  proce¬ 
dures.  Subsequent  discussion  of  Q-values  here  refer  to  those  defined  by  (2.9). 

Figure  4.7  summarizes  the  J-Q  description  of  stress  triaxiality  at  different  points  along 
the  crack  firont.  Both  the  viniaxial  and  biaxial  cases  are  taken  to  the  same  load  level,  i.e., 
ct*/cto=1.1-  Q  is  evaluated  at  =  2ahead  of  the  crack  front.  Under  uniaxial  load¬ 

ing,  Q  values  for<^  >  45”  saturate  at -0.8  for  large-scale  yielding.  Near  the  free  surface,  0=0, 
stress  triaxiality  is  reduced  to  a  level  approaching  the  yield  stress  even  at  relatively  low 
loads  (Q-*-2,0).  Biaxial  loading  promotes  essentially  uniform  stress  triaxiality,  Q=-0.7, 
overmuch  of  the  crack  front.  However,  at  the  0=2.4”  and  17”  locations,  the  influence  of  biax¬ 
ial  loading  is  very  pronounced.  Final  Q-values  for  these  two  locations  reveal  an  increase 
in  stress  triaxiality  on  the  order  of  the  yield  stress  relative  to  the  uniaxial  loading  response. 


Q  Q 

Figure  4.7  Evolution  of  stress  triaxiality  in  surface-cracked  plate  with  increasing  local 
deformation,  as  measured  hy  Jiocallf^o,  along  radial  lines  emanating  from  dif¬ 
ferent  points  on  the  crack  front.  Q  is  evaluated  at  r/(J/oco/(oo)=2.  Both  geome¬ 
tries  are  loaded  to  the  same  level  o»/ao=l.l. 

For  crack  front  locations  0-*O,  the  mechanically  imposed  biaxial  stress,  o^,  corre¬ 
sponds  to  a  positive  T-stress  (a  stress  parallel  to.X^,  see  Fig.  4.1).  Under  SSY  conditions 
in  the  surface-cracked  plate,  the  T-stress  elevates  Q  slightly  above  zero  in  accord  with  the 
discussion  in  Section  2.  Under  LSY  conditions,  the  T-stress  brings  about  a  higher  level  of 
stress  triaxiality  near  the  free  surface;  however,  the  Q-values  are  still  negative  indicating 
a  loss  of  stress  triaxiality  relative  to  the  high  constraint,  reference  condition  of  plane-strain 
SSYt*o-  Thus,  at  an  identical  value  of  applied-J  in  uniaxial  and  biaxial  loading,  the  the 
crack  front  region  with  maximum  opening  mode  stress  occurs  near  0  *=  17”  for  the  biaxial 
loading.  However,  the  magnitude  of  remote  loading  required  to  generate  the  equivalent  J- 
values  is  larger  for  the  biaxial  case  (see  Fig.  4.3  and  4.4;  note  the  overall  larger  J-vedues 
for  viniaxial  loading). 


28 


4.3  Matching  Structural  and  Test  Specimen  Constraint 


Consider  the  fracture  assessment  of  a  structural  configuration  which  is  modeled  adequate¬ 
ly  by  the  surface  cracked  plate  subjected  to  uniaxial  or  biaxial  loading.  The  J-Q  and  micro¬ 
mechanics  concepts  provide  quantitative  frameworks  to  select  a  corresponding  laboratory 
test  specimen,  an  SE(B)  for  example,  which  produces  the  same  crack  front  constraint  as  the 
structural  configuration.  The  cleavage  fracture  toughness,  Jc,  measured  with  such  a  speci¬ 
men  should  then  be  employed  in  fracture  assessments  of  the  structure.  These  two  ap¬ 
proaches  are  illustrated  in  Figs.  4.8  and  4.9. 

Figure  4.8  compares  the  computed  J-Q  drivingforce  curves  for  SE(B)  specimens  having 
a  range  of  a/ W  ratios  with  the  driving  force  curves  at  0=90°  (uniaxial  loading)  and  0=17° 
(biaxial  loading)  for  the  siuface  cracked  plate.  An  SE(B)  specimen  with  a/WsO.OS  best 
matches  the  evolution  of  stress  triaxiality  for  the  uniaxially  loaded  plate  while  an 
afW «0.20  best  matches  triaxiality  for  the  biaxial  loading.  By  using  the  J-Q  description  of 
crack  front  stresses  as  input  to  the  micromechanics  model  (as  described  in  Section  3.4),  the 
effects  of  constraint  on  cleavage  fracture  toughness  for  the  surface  cracked  plate  and  SE(B) 
specimens  of  selected  ratios  are  predicted  as  shown  in  Fig.  4.9.  An  SE(B)  specimen  with 
a/W=0.05  very  closely  matches  the  uniaxial  loading  curve  for  0=90°  while  an  a/W*0.20 
SE(B)  specimen  closely  matches  the  biaxial  loading  curve  for  0=17°. 

Table  4.1  shu  ws  predicted  values  of  the  remote  stress  at  cleavage  fracture  determined 
by  the  mic.omechanics  approach.  Fracture  occurs  when  the  applied  Jo  reaches  a  critical 
value  at  the  crack  front  position  exhibiting  the  most  stress  triaxiality,  i.e.,  the  locations  at 
which  Jiocal  generates  the  maximum  Jo  values.  As  indicated  by  Figs.  4.7  and  4.9,  critical 
locations  on  the  crack  front  are  0=17°  for  biaxial  loading  and  0=90°  for  uniaxial  loading. 
In  Table  4.1,  three  representative  values  of  Jo/(aao)  at  fracture  are  considered.  Correspond¬ 
ing  values  of  Jiocal  for  the  biaxial  and  uniaxial  loadings  are  fovmd  from  Fig.  4.9  wi^  the 
applied  loads  to  produce  these  J-values  given  by  Figs.  4.3  and  4.4.  These  computations  sug¬ 
gest  that,  despite  the  higher  stress  triaxiality  of  the  biaxial  loading  at  0= 1 7°,  slightly  larger 
fracture  loads  are  predicted  for  the  biaxial  loading.  For  this  combination  of  material  proper¬ 
ties  and  crack-specimen  geometry,  the  effects  of  higher  stress  triaxiality  at0=17°  are  offset 
by  the  lower  applied  J-values  at  this  location  on  the  crack  front  relative  to  the  0=90°  loca¬ 
tion. 

The  potential  advantage  offered  by  the  micromechanics  approach  becomes  clear  from 
Fig.  4.9.  It  is  not  necessary  to  determine  which  laboratory  specimen  matches  the  structural 
constraint;  rather,  any  afW  ratio  SE(B)  can  be  tested  to  measure  the  size  independent  frac¬ 
ture  toughness,  Jo,  from  which  the  structural  toughness,  Jc,  for  each  loading  case  (uniaxial 
and  biaxial)  is  predicted  from  the  corresponding  structural  response  curves  shown  in  Fig. 
4.9.  When  deep  notch  SE(B)  data  is  already  available,  from  material  qualification  tests  for 
example,  no  additional  shallow-crack  testing  is  needed  to  apply  the  micromechanics  ap¬ 
proach.  Such  applications  of  the  micromechanics  model  imply  that  the  same  fracture  mode 
(cleavage)  occurs  in  both  the  laboratory  specimen  and  the  structural  configuration.  The  cvir- 
rent  model  cannot  predict  the  effects  of  specimen  geometry  and  loading  mode  on  the  frac¬ 
ture  toughness  (J/d  which  characterizes  the  initiation  of  stable,  ductile  tearing.  However, 
the  model  does  predict  when  fracture  by  cleavage  becomes  highly  unlikely.  Consider  the 
response  for  the  o/W=0.5  SE(B)  specimen  shown  in  Fig.  4.9;  if  the  Jo  measured  with  this 
specimen  is  sufficiently  large,  the  driving  force  curve  for  a  shallow  notch  SE(B)  specimen 
or  for  the  surface-cracked  plate  never  attain  such  a  large  value  of  Jo.  The  model  predicts 
that  cleavage,  without  prior  tearing,  does  not  occur,  i.e.,  the  interaction  of  crack-tip  plastic 
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Q 

Figure  4.8  Application  of  J~Q  methodology  to  select  a/W  ratio  in  SE(B)  specimen  for  test¬ 

ing  that  matches  stress  triaxiality  in  structural  configuration  represented  by 


Figure  4.9  Application  of  micromechanics  methodology  to  select  a/W  ratio  in  SE(B)  spec¬ 

imen  for  testing  that  matches  effective  driving  force  for  cleavage  fracture  in 
structural  configuration  represented  by  a  surface  cracked  plate  loading  in 
uniaxial  or  biaxial  tension. 
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zones  and  nearby  free  surfaces  prevents  near-tip  stresses  from  achieving  the  critical  levels 
needed  to  trigger  cleavage. 

The  SE(B)  responses  employed  in  this  discussion  are  obtained  from  2-D,  plane-strain 
computations.  In  such  cases,  constraint  matching  with  a  structural  configuration  is  accom¬ 
plished  by  var3dng  the  absolute  specimen  size  and/or  the  alW  ratio.  For  the  same  a/W  ratio, 
large  specimens  increase  constraint  at  a  given  J-value  relative  to  small  specimens.  Simi¬ 
larly,  for  a  fixed  specimen  size,  large  afV/  ratios  increase  constraint  at  a  given  J-value  rela¬ 
tive  to  small  a/W  ratios.  Different  thicknesses  provide  yet  another  means  to  vary  constraint 
in  test  specimens,  although  experimental  and  computational  evidence  suggest  the  thick¬ 
ness  effect  is  much  less  significant  than  absolute  size  or  a/W  effects  when  specimens  of  usu¬ 
al  proportions  are  employed  (S  >  W/2). 


Table  4.1:  Predicted  Loads  at  Fracture  for  Surface-Cracked  Plate 

A 

CUJq 

Biaxial  Loading  (^=17°) 

Uniaxial  Loading  (0=90°) 

local 

CUJq 

m 

^ local 

CT" 

(^0 

0.00333 

0.00549 

0.69 

0.89 

0.00960 

1.20 

0.86 

0.00666 

0.01453 

1.82 

1.06 

0.02650 

3.31 

1.01 

0.00833 

0.02300 

2.88 

1.14 

0.03682 

4.60 

1.05 
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5.  CONCLUSIONS 


Our  investigations  have  shown  that  two-parameters,  J  and  Q,  suffice  to  characterize  the 
full  range  of  near-tip  environments  at  the  onset  of  fracture.  J  sets  the  size  scale  of  the  zone 
of  high  stresses  and  large  deformations  while  Q  scales  the  near-tip  stress  level  relative  to 
a  high  triaxiality  reference  stress  state.  The  structure  of  the  J-Q  fields  has  been  estab¬ 
lished  by  higher-order  asypmtotic  analysis  and  full-field  numerical  calculations  within 
the  context  of  the  modified  boundary  layer  formulation.  Detailed  analyses  of  finite-width, 
crack  bodies  show  that  the  J-Q  fields  dominate  over  physically  significant  size  scales,  i.e. 
they  represent  the  environment  in  which  the  ductile  and  brittle  mechanisms  are  operative. 
Therefore,  the  J-Q  fields  furnish  the  theoretical  basis  to  address  onset  of  clevage  fracture, 
the  initiation  of  ductile  tearing,  as  well  as  the  competiticn  between  cleavage  fracture  and 
ductile  tearing.  Indeed,  the  J-Q  theory  can  provide  a  framework  which  allows  the  cleavage 
and  ductile  toughness  loci  to  be  measured  and  utilized  in  engineering  applications. 

Constraint  effects  on  cleavage  fracture  have  been  the  subject  of  a  number  of  recent  stud¬ 
ies.  The  J-Q  theory  together  with  a  micromechanical  model  for  cleavage  predicts  that 
cleavage  fracture  toughness  depends  sensitively  on  near-tip  stress  triaxiality.  The  cleav¬ 
age  toughness  locus  has  been  measured,  for  example  by  Sumpter  and  Forbes  [10]  for  a  mild 
steel  and  by  Kirk,  et  al.  [42]  for  A515  steel.  The  toughness  data  do  show  a  strong 
dependence  on  Q.  Tbughness  elevations  of  about  5  or  so  have  been  measured  in  low 
constraint  crack  geometries.  Constraint  is  also  expected  to  exert  an  influence  on  the 
initiation  of  ductile  tearing;  however,  mechanistic  studies  of  ductile  tearing  and  the  limited 
experimental  data  suggest  that  ductile  initiation  toughness  depends  less  strongly  on  stress 
triaxiality.  Systematic  experimental  studies  are  required  to  quantify  constraint  effects  on 
the  initiation  of  ductile  tearing. 

The  experimental  determination  of  a  toughness  locus  can  become  very  costly,  requiring 
considerable  material  and  testing  time,  especially  if  toughness  data  are  required  for  a 
number  of  temperatures.  An  alternative  approach  for  cleavage  fracture  appears  feasible. 
The  constraint  correction  procedure  advocated  here  uses  a  limited  experimental  database 
to  predict  cleavage  toughness  over  a  broad  range  of  stress  triaxiahty.  The  procedure  has 
been  applied  to  several  series  experiments  and  the  results  are  very  encouraging.  The 
procedure  and  its  theoretical  basis  are  discussed  in  Section  3.  Though  it  has  not  been 
discussed,  a  similar  procedure  can  be  developed  (in  principle)  for  the  initiation  of  ductile 
tearing.  As  is  the  case  for  cleavage  fracture  toughness,  the  procedure  can  be  used  in 
conjvmction  with  the  J-Q  fields  to  determine  geometry-  and  load-dependent  ductile 
fracture  toughness  data. 

Procedvures  such  as  the  one  in  Section  3  hold  promise  for  predicting  toughness  loci  for 
cleavage  and  ductile  fracture  and  can  facilitate  engineering  applications  of  the  J-Q  ap¬ 
proach.  However,  the  incorporation  of  micromechanics  failure  criterion  in  a  fracture  me¬ 
chanics  methodology  is  not  without  its  limitations.  It  is  essentiEd  that  the  failure  criterion 
be  a  sufficiently  realistic  model  of  the  actual  fracture  process.  In  applications  where 
cleavage  and  ductile  failure  modes  are  operative,  competing  failure  processes  also  must  be 
considered.  Unfortunately,  the  initiation  of  ductile  tearing  is  adso  a  process  involving 
several  competing  mechanisms,  such  as  void  formation  vs.  void  growth  and  coalescence, 
void  sheet  formation  and  shear  localization,  and  as  yet  there  is  no  general  agreement  as 
to  the  essential  features  of  a  realistic  model.  Micromechanicsd  models  for  ductile  tearing 
are  necessarily  more  complex,  involving  more  metallurgical  properties  for  a  material,  than 
the  models  for  cleavage  fracture.  Considering  the  enormous  economic  payoffs,  however, 
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greater  efforts  are  warranted  to  establish  realistic,  robust  micromechanical  models  for 
ductile  fracture. 
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i3.  AB5TKAC»  i^^aximu'Ti  200  oras, 

This  report  explores  the  fundamental  concepts  of  the  J-Q  description  of  crack-tip 
fields/  the  fracture  toughness  locus  and  microroechanics  approaches  to  predict  the 
variability  of  macroscopic  fracture  toughness  with  constraint  under  eleistic-plastic 
conditions.  While  these  concepts  derived  from  plane-strain  considerations,  initial 
applications  in  fully  3Br>geoenetries  are  very  promising.  Computational  results  are 
presented  for  a  surface  cracked  plate  containing  a  6:1  sani -elliptical ,  a=t|4  flaw 
subjected  to  remote  uniaxial  and  biaxial  tension.  Crack-tip  stress  fields  consis¬ 
tent  with  the  J-Q  theory  are  demonstrated  to  exit  at  each  location  along  the  crack 
front.  The  microroechanics  model  arploys  the  J-Q  description  of  crack-front  stress 
to  interpret  fracture  toughness  values  measured  on  le±)oratory  specimens  for  fract¬ 
ure  assessment  of  the  surface  cracked  plate.  The  computational  results  suggest 
only  a  minor  effect  of  the  biaxial  loading  on  the  crack  tip  stress  fields  and,  con¬ 
sequently,  on  the  propensity  for  fracture  relative  to  the  uniaxial  loading. 
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